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^ . Abstract 

■^j- ■ The light curve of a type la supernova decays at a rate set by the /3-decay lifetimes 

of the 56 Ni and 56 Co produced in the explosion. This makes such a light curve sensitive 
to the value of the Fermi constant Gf at the time of the supernova. Using data from 
the CfA Supernova Archive, we measure the dependence of the light curve decay rate on 



redshift and place a bound on the time variation of Gf of 



< 1Q- 9 yr 



1 Introduction 



Finding new fundamental physics beyond the standard model and general relativity has 
long been a topic of great interest. One way to identify definitive evidence of new physics 
is to look for exotic phenomena that simply cannot occur in the standard physical theories. 
One such exotic possibility that has provoked significant interest is having time-dependent 

fundamental constants [UI21E1HI51E1IZ1I8]. 

The possibility of time-varying fundamental constants opens up many new lines of 
inquiry. The most commonly studied theory with varying constants is one with a time- 
dependent fine structure constant a. A quantum approach to this subject can be found 
in [9], while other approaches treat the problem rather phenomenalistically (e.g. [10J). 
Another possibility that has been fairly well studied is a change in the quantum chromo- 
dynamics scale A.qcd and therefore the electron-proton mass ratio. 

If some of the fundamental constants in nature are allowed to change over time, the 
effects of these changes might be detected. Searches for such temporal variations entail 
measurements of the same experimental observable at different points in time, to see 
if there is any change. However, since the changes involved should be tiny, laboratory 
experiments searching for them demand very high precision. Less precision is required if 
large time scales are involved; a very slow time variation can be observed if the times used 
to detect it are long enough. Thus, cosmological observations, which give us a view of 
how physics operated millions or billions of years in the past, seem to be good candidates 
to test for time variations in fundamental constants. 

The purpose of this paper is to study a possible change in the Fermi constant Gp. 
This quantity is defined to be Gp = with g the standard model's SU(2)l gauge 

coupling constant and Mw the mass of the W boson. Since the vector boson mass is 
Mw — hgv, Gp is really a measure of v, the Higgs field vacuum expectation value. The 
masses of fundamental fermions depend on v and the values of the fermion-Higgs Yukawa 
couplings Xf, with rrif = XfV. So a changing value of Gp would also affect the fermion 
masses unless the couplings Xf also changed to compensate. Since there is no evidence 
of any time variation in m e in atomic spectra, we shall assume that such cancellation 
between the time variation of v and Xf occurs, even though that is not a particularly 
natural scenario. 

Gp also governs the strength of charged current weak interactions at low energies. It 
can be measured with relatively high precision by studying the lifetimes of muons, and 
its accepted value is Gp = (1.6637 ± 0.0001) x 10~ 5 GeV" 2 [EEJ [Lj. However, laboratory 
measurements of decay rates are not precise enough to offer a particularly promising way 
of measuring Gp. On the other hand, the light coming from supernovae explosions can 
provide a useful way to study weak interactions that occurred at much earlier times. 

In Section we shall discuss the physics of supernovae and how their light curves are 
related to the Fermi coupling constant Gp. Section [3] covers the details of the calculational 
methodology used to relate the supernova light curve data to Gp. In Section HI we apply 
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this methodology to data from the Harvard-Smithsonian Center for Astrophysics (CfA) 
Supernova Archive [HJ US] and state our results. Section [5] presents conclusions and 
outlook. 

2 Supernova Structure 

Supernovae are catastrophic variable stars that, having previously produced energy through 
steady-state fusion processes, explode and release huge additional amounts of energy. The 
amount of energy released is so vast that even cosmologically distant supernovae can be 
easily detected in the sky. Normal stars avoid gravitational collapse because nuclear fusion 
processes are constantly generating energy in their cores, where the high temperature and 
pressure enable light nuclei to overcome their electrostatic repulsion and fuse. A large red 
giant star may burn a significant portion of its mass into isotopes of iron, beyond which 
further fusion is no longer exothermic. However, smaller stars, without thick external lay- 
ers to compress the core to the high temperatures and pressures needed to produce iron, 
may cease burning when the core is mostly carbon and oxygen, without having extracted 
all of this fuel's available energy. In either case, the core of the star settles down to become 
a white dwarf, supported against gravitational collapse by the degeneracy pressure of its 
electrons. 

If the mass of a white dwarf star is larger than approximately 1.4M Q (the Chan- 
drasekhar mass Mch), the electron degeneracy pressure can support the star against 
gravitational collapse. However, if the mass of the star grows larger than this limit, the 
electron pressure can no longer prevent a collapse. Such a collapse triggers the release of 
a tremendous amount of energy and the star becomes a supernova. Iron white dwarfs at 
the cores of massive stars collapse to form even more compact objects — neutron stars or 
black holes. However, of interest here are the type la supernovae produced from carbon- 
oxygen white dwarfs. These stars may accrete mass from (or possibly merge with) their 
partners in binary systems, until the white dwarf's mass becomes sufficient to trigger a 
supernova. When these stars exceed Mch, they begin to collapse. However, the collapse 
phase is short-lived; as the star contracts, there is an increase in temperature and density 
in the core of the star, which allows the ignition of carbon fusion. Within a short period of 
time, a huge amount of fusion energy is released, blowing the star apart and increasing its 
luminosity to enormous values. (The total luminosity may be that of millions of ordinary 
suns) . 

In this paper, we shall be analyzing the light curves of type la supernovae. The near 
equality of the masses of the white dwarf progenitors gives them all similar energy outputs 
and makes them extremely important standard candles in cosmology. Type la supernovae 
all seem to follow a characteristic light curve. (Some corrections to the light curves are 
required in order for the population of type la supernovae to act as genuinely standard 
candles, but these are not important for our present purposes.) The relative uniformity 
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of this curve structure for different stars also makes them a very useful tool for studying 
certain physical properties. 

The light curve is a plot of the apparent magnitude m = —2.512 log -j- of the source 
versus time; / is the radiation flux from the source, and /q is a reference value of the 
flux defining magnitude. Using different kinds of filters, which cover different parts of 
the spectrum, several different magnitudes may be measured for a single source. There 
is a separate light curve for each band, but the curves for different colors show similar 
behavior, which is again consistent across different sources. 

At early times, the brightness of a supernova increases rapidly until it reaches a max- 
imum, which is followed by two regions in which m decreases linearly (corresponding to 
exponential decays in the flux). Since type la supernovae are ultimately powered by un- 
controlled nuclear fusion, they produce large quantities of the heavy elements, especially 
those elements with binding energies close to the maximum value. Specifically, each su- 
pernova produces a great deal of the double magic nucleus 56 Ni. This nucleus decays by 
electron capture (EC) to 56 Co with a lifetime of 6.077 d, and this decay provides much 
of the energy released by the supernova during its earlier stages. The first linear region 
along the light curve corresponds to the period during which 56 Ni decay is the dominant 
source of heating. The daughter nucleus 56 Co also decays by EC to 56 Fe, with a lifetime of 
77.27 d. At later times (in the second, more slowly decaying linear region), it is this decay 
that primarily fuels the supernova's luminosity. This region, whose extent is determined 
by the decay rates of 56 Ni and 56 Co, can be found between approximately 35 and 100 
days after the maximum magnitude is reached. During this time, some months after a 
supernova reaches its peak, the light curve approximately tracks the quantity of 56 Co still 
present. The falloff in the luminosity can therefore be used to measure the lifetime of this 
nuclide at the time the supernova occurred. This gives direct information about the value 
of several standard model parameters at the time of the explosion's occurrence, including 
Gf- This was first suggested in [13], when much less data was available. 

An EC decay rate depends on fundamental constants in the following fashion. Being a 
tree- level weak decay, the rate is naturally proportional to G 2 F . However, it also depends 
on the probability that an s-state electron will be found close enough to the nucleus to 
be taken up by a proton. The nuclear volume is proportional to Aq^ d . The square of the 
electron wave function in the vicinity of the nucleus is proportional to 3 , where ao = 
is the Bohr radius. The inverse V of the mean life of a nucleus that is unstable against EC 



(The rate also depends on the Q-value of the decay, which is further dependent on a 
and A ™* D , in more complicated fashion.) However, changes in a and A ^ D oc ^ are 
much easier to measure than would be changes in Gp. ct and ^ affect atomic energy 
levels directly, so their time derivatives can be constrained by comparing spectra taken 
at different times. Cosmological and laboratory measurements have constrained these 
dimensionless quantities to vary only by parts in 10 16 or less per year. 



is therefore proportional to the combination of fundamental 
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A supernova's light curve and total luminosity are also affected by the value of the 
gravitational constant GV- The amount of 56 Ni produced in the explosion is proportional 
to the total mass of the collapsing star, M Ch oc G^ 3 ^ 2 [TBI E] • Observations of supernovae 



have been used to place bounds on the variation of G^ at the 



< 10 -11 yr" 1 level. 



We choose here to look at a model that neglects variations in Gn, although an analysis 
that considered simultaneously varying Gn and Gp would also be interesting. We shall 
also neglect any possible time dependences of a and , in light of the strong bounds 

on such time dependences. This means there is a simple relationship between T and Gp, 

^ = --. (1) 
Gf 2 r K J 

Bounds on changes in the 56 Co decay rate may thus be interpreted as bounds on a change 
in the Fermi constant. 



3 Method 

In the light curve of a type la supernova, the slope of the second linear region is related 
to the decay rate of 56 Co. The main goal of this paper is to compare this slope across 
different supernovae of different ages. The slope for each one is related to the Fermi 
constant at a particular time. We will make this comparison using data from the CfA 
Supernova Archive. 

Ideally, one could simply measure the slopes of the supernova light curves in their 
56 Co-dominated regions and then determine whether the slopes are correlated with the 
ages of the events. However, several corrections must be made before this can be done. 

The light curves do not give a pure measurement of the 56 Co lifetime; what is observed 
is the luminosity of the supernova remnant. The EC decay is the object's main energy 
source, but the decay energy is moderated through the complicated interactions of the 
expanding gas, so the luminosity does not track perfectly with the 56 Co decay profile. 
The details of the remnant dynamics differ a bit from supernova to supernova. It is not 
possible to correct for this lack of complete uniformity in the structures, and this will 
necessarily be a source of random error. 

Some of the corrections may be handled more simply, however. This is the case for 
the correction arising from the well-known Doppler shift associated with the expansion of 
the Universe. Since the inverse of the mean life, T, measures how fast 56 Co is decaying 
(and therefore, the frequency at which EC processes are occurring), the observed value 
of T will be subject to the same cosmological redshift as any other frequency. In terms 
of the red shift z, the measured value T on Earth and the proper decay rate T at the 
exploding star are related by 

r„ = (i + *)r. (2) 
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The physical effect of the redshift is to slow down the decay of the light curve. 

We shall denote by m the slope of the second linear region of each particular light 
curve, as observed on the Earth. This slope and the 56 Co decay rate should be directly 
proportional, so we have 

r„ = (1 + z)Y = A(l + z)m = Am , (3) 

where mo is the proper decay rate of the supernova light curve, and A is a proportionality 
constant related to the complex physics of the explosion. The quantity mo is proportional 
to the square of the Fermi coupling constant at the time t Q that the explosion occurred, 
m oc T oc G F (t ). If Gp has changed over time, we would expect to see a change in 
the decay rates m observed for different supernovae. Comparing the light curves for 
supernovae of different ages gives a good way to measure a change in the Fermi coupling 
constant. 

The light from a particular supernova is emitted and collected at the times to and t, 
respectively. Denoting the scale factor of the Universe by a(t) and the Hubble constant 
by H, we have the standard relations 

a(t ) 

H = ». (5) 
a 

The factor H converts derivatives with respect to z into ones with respect to time. Based 
solely on the relations (jlj) and (J5J), one may be tempted to conclude that j£ = —j^^j^- 
(where (rho) is the average value of rh as computed from the different sources). However, 
a new complication arises because of the differences in the behaviors of the m observed 
in different wavelength bands, and this relation for ^ is subject to further corrections. 
Looking at the light from different color regions, we see slightly different slopes associated 
with the different bands. At first glance this might not seem to be a problem, since the 
value of the proportionality constant A defined in eq. (j3]) could just be redefined for each 
band. The value ^ does not depend on A, so it would not affect the final result. 

Nevertheless, the fact that A does not remain constant across the various bands in- 
troduces a new correction. To understand why a further correction is necessary, let us 
consider light that is detected on Earth using a red filter. This light has been redshifted 
by the expansion of the Universe. If the slopes for the different bands were the same 
this would not be a problem, but the red light that is detected actually had a shorter 
wavelength when it was emitted during the explosion. All the red light that is detected on 
Earth was some amount bluer at the time of its emission; and we should ideally be compar- 
ing light from different sources that had the same wavelength at the time of emission — not 
at the time of observation. If the light curve steepens with increasing wavevector k, the 
redshift would transform this into a steepening (within a given wavelength band) as a 
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function of z, because the older light originally had a higher frequency. What is required 
is similar to a light curve K correction [18]. However, rather than using tabulated K 
corrections, we shall, for this low redshift sample of supernovae, calculate our correction 
directly from the data set. This is justified because the technique minimizes our reliance 
on the results of unrelated fits and because we shall find that the correction it produces 
is relatively small. 

To account for the redshift 's effect, we suppose that the dependence of ttlq on k is 
weak, and make a linear approximation 

(1 + z)m x = (m x )o + ~J]r( k o ~ k *) = (™*)o + ~^^ zkx ">- ( 6 ) 

% is the observed slope in the band X with mean wavenumber kx- Once corrected 
by (1 + z), nix gives the proper decay rate m , but for a band centered at a higher 
wavenumber k = (1 + z)k X - The z dependence of k generates a spurious dependence on 
z in the observational data, which we must subtract away. 

Unlike the simple (1 + z) factor relating m to mo (which is an exact correction), there 
is no particular reason for this effect to be linear in z. The wavelength dependence of the 
proportionality constant A is a consequence of the complicated dynamics of the supernova 
remnants; its detailed form is not known. If is itself a function of z, there will be 
0(z 2 ) corrections to (EJ). At the same order, it would also be necessary to correct the 
average ((rhx)o) to account for the fact that we are actually averaging decay rates in 
different color regions. We are neglecting these possibilities (which would correspond to 
corrections that are higher order in z in our final formula for ^), because the redshifts 
in the sample we shall study range only from 3.9 x 1CT 3 to 2.4 x 10 -2 . However, these 
corrections can be computed, and working with sources at higher z would require us to 
have a better knowledge of how rho depends on k. Yet there is a limit to the precision 
that can be derived from more detailed models, since there are data available for only a 
handful of wavelength bands (five in this case), and data points at multiple frequencies 
would be required to map out the detailed dependence of A on k. 

In principle, we could try to correct the values of rho for each supernova individually, by 
studying the dependence of its light curves on k. We could then use the corrected values 
obtained separated for each source to compute j^^r- However, this would require 
additional fitting for each supernova, which may introduce additional errors and which 
cannot be practically implemented for many of the sources because of the limited data 
available. We shall therefore make only a single calculation of ^ , using the combined 
data from all the sources. 



Using eq. ([6]) we find our final formula for ^ 



r h 



r (mo) 



dfn,Q 2n f d(rho) 



(7) 



dz A \ dk 

including the linear correction to eliminate the spurious z dependence caused by a nonzero 

drag 
dk 



^jf-. Evaluating this requires multiple fits. The order of the fitting procedures is as follows. 
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m-x 



♦ m v 

x m B -2 
m v -4 
a m R -6 
□ mi -8 





Figure 1: Light curves for SN 1998aq. The magnitudes for the different bands are mjj, 
tub — 2, my — 4, rriR — 6, and mi — 8. The black lines represent linear fits between 50 
and 130 days after the peak. 



First, rho is determined for each supernova in each band separately. These are then fitted 
as a linear function of the redshift z, giving both an average (rh ) and a slope ^p 3 - for 
each band. From the fit of the five average values (too) versus k, the parameter d ^ can 
be extracted. Then j£ can be evaluated according to eq. (JBj); this evaluation can be done 

separately in each of the five bands, although the single fitted value for the quantity 
is common to all of them. 

The overall negative sign in eq. ([7]) merits some explanation. Because the light curves 
are decaying in the region of interest, tjiq is always negative. If tjiq is more negative at 
greater values of z < and j^^r > 0), the decay rate was faster at earlier times. 
This means r and thus Gp are decreasing as time passes, leading to the negative sign in 
eq. 0. 



4 Data Analysis 

Using the method we have discussed, we can now place a bound on the time variation rate 
of Gp using data from the CfA Supernova Archive. For purposes of illustration, we shall 
initially focus on one particular supernova, SN 1998aq [19] . for which z = 3.95 x 10 -3 . Its 
light curves are shown in fig. [TJ plotting the supernova's magnitude as a function of time t 
in days. The linear behavior arising from the 56 Co decay can be seen in a region covering 
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Band 


m (lO^d- 1 ) 


U 


-214 ±7 


B 


-138 ± 2 


V 


-261 ±4 


R 


-321 ±4 


I 


-415 ±6 



Table 1: Fitted values of the slope fn for SN1998aq, in units of 10 4 d 1 . 

approximately 50-130 days after the maximum magnitude was recorded. Although the 
data has small errors related to the measurements, these have been omitted for clarity of 
presentation. To find the slopes in the different bands, we used a weighted linear fit in the 
region of interest, minimizing the x 2 per degree of freedom. The results of these fits are 
also shown in fig. [T] and the fitted slopes are listed in table [TJ In order to transform these 
results into the reference frame of supernova, we just use the relation rho = (1 + z)m. 

We have applied the same procedure to multiple sources. Not all sources in the 
CfA Supernova Archive had sufficient data to be fitted reliably, and not all bands were 
represented for all sources. In table (2j we summarize the results found for the different 
supernovae we analyzed. The slopes m are in the observer reference frame on Earth (no 
1 + z factor included). There was one further complication associated with fitting some 
of the light curves. For some sources, there is limited data at early times, and it was not 
possible to resolve the date at which the light curve achieved its absolute maximum. In 
many cases, we could identify the proper linear region by inspection, but if there was any 
doubt whether a given data point should be included, the fitted region was chosen so as 
to produce the best linear fit. 

The data available for some supernovae and some wavelength bands were not of suffi- 
cient quality to be useful. We omitted those sources for which the standard error in the 
fit of fn was larger than 0.015 d -1 and those with fewer than three data points in the re- 
gion of interest. The former was indicative of substantial deviations from linear behavior 
(that is, nonexponential decay in the luminosity), and the later simply represented too 
few points to be useful. 

The z dependence of fn for each wavelength band is shown in fig. [2j It is evident that 
^p- is negative in each of the five distinct wavelength regimes. The values of the slopes, 
as well as some other relevant quantities, are collected in table [3J 

The last fit required was for the wavelength dependence of (fn Q ). The result was 
2ir d<i ™^ = (9.3 ± 4.8) nm/d. The correction associated with this factor will make only a 
small contribution to the final value of This correction, along with the other factors 

necessary to produce the final computed value of j^, are given in table HI The fourth 
column of this table collects all the corrections for each band, giving the full z-dependence 
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(m )u (1(T 2 d" 1 ) (m ) B (10" d-') 




(mo), (1(T 2 d"') 




t , , , , z 

0.005 0.01 0.015 0.02 

Figure 2: Plots of m versus z for the U, B, V, R, and I bands (denoted by colors as 
in fig. [1]). In each plot, the dashed horizontal line marks the mean (mo); the substantial 
difference in this mean for different wavelength regions indicates a nonzero . 
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Source 


10 5 z 


mv (10 _4 d _1 ) 


rh B (10 _4 d _1 ) 


my (10~ 4 d _1 ) 


m R (10 _4 d _1 ) 


mi (10 _4 d~ 1 ) 


1 QQ^Kr] 

lyyoDQ 


i fin 




1 71 _i_ 7 


971 _1_ Q 






1998aq 


395 


—214 ± 8 


— 138 ± 2 


—261 ± 4 


— 321 ± 4 


—415 ± 6 


1994ae 


427 




-145 ±2 


-224 ±4 


-286 ± 7 


-359 ± 8 


1995al 


514 




-152 ±8 


-246 ±4 


-311 ± 3 


-366 ± 1 


1995D 


654 




-131 ± 5 


-248 ± 6 


-333 ±35 


-384 ± 24 


1999gh 


767 


-182 ± 20 


-175 ±5 


-251 ± 5 


-322 ± 12 


-299 ± 33 


2000cx 


793 


-226 ±4 


-187 ±2 


-257 ± 8 


-312 ±8 


-347 ± 19 


1999ac 


949 


-240 ± 6 


-175 ± 14 


-231 ± 11 


-313 ± 13 


-317 ±40 


1998bp 


1042 




-187 ± 13 






-289 ± 42 


1998es 


1056 




-187 ± 16 


-314 ± 15 


-376 ± 7 




1995E 


1156 






-253 ± 20 


-321 ± 17 


-441 ± 23 


1999dq 


1432 


-254 ± 1 


-163 ± 15 


-250 ± 17 


-332 ± 18 


-387 ± 33 


1999aa 


1443 


-210 ±2 


-130 ±8 


-248 ± 7 


-316 ±6 


-414 ± 7 


1993ae 


1904 






-231 ± 41 


-316 ± 24 


-428 ± 123 


1994M 


2296 




-149 ±7 


-294 ± 18 


-399 ± 2 


-454 ± 1 


1995ak 


2300 




-243 ± 20 


-310 ± 12 


-418 ± 25 


-489 ± 57 



Table 2: Fitted values of m, in units of 10 4 d 1 . Sources and bands for which the data 
were insufficient are marked "— " . 



Band 


A (nm) 


^ (io-M- 1 ) 


(rho) (lO^d" 1 ) 


U 


350.0 


-23 ± 27 


-22 ±3 


B 


438.0 


-19 ± 12 


-17 ±3 


V 


546.5 


-28 ± 10 


-26 ±3 


R 


647.0 


-46 ± 11 


-34 ±4 


I 


786.5 


-52 ± 19 


-39 ±7 



Table 3: Fitted means and ^-dependences of rho for the five wavelength bands. A is the 
central wavelength for each band. 

of the 56 Co decay rate T ; only a factor of the Hubble constant H is needed to convert 
this into a time dependence. 

Using the values from table HI we may compute the weighted average of ^ for all the 

frequency bands. The result is = — 12.7 ± 7.4, so 

^ = (-4.8 ± 2.8) x 10- 10 yr- 1 . (8) 

Gr p 

where we have used H = (74.2 ±3.6) km/s/Mpc [20], the uncertainty in which introduces 
a small additional error. 
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m v 




t(d) 



20 40 60 80 



Figure 3: V-band light curves of SN 1995ak and SN 1998aq, showing linear fits in the 
56 Co decay region. The maximum intensity occurs at t = 0, where my is normalized to 
— 12.5 and —12, respectively. 



Band 


1 drriQ 
(mo) dz 


k d{mo) 
(mo) dk 


1 f o 
H To 


U 


-10.5 ± 12.1 


-1.2 ±0.6 


— 1 1 .7 ± 12.1 


B 


-11.5 ±7.7 


-1.3 ±0.7 


— 12.8 ± 7.8 


V 


-10.6 ±4.0 


-0.7 ±0.3 


-11.2 ±4.0 


R 


-13.6 ±3.8 


-0.4 ±0.2 


-14.1 ±3.8 


I 


-13.4 ±5.5 


-0.3 ±0.2 


-13.7±5.5 



Table 4: Contributions to using supernovae in the domain 3.9 x 10 3 < z < 2.4 x 10 2 . 

5 Conclusions 

Our main result, eq. OH]), shows that the Fermi constant Gp may be decreasing with time; 
however, the evidence for a change is at less than the 2a level. The error is too large to 
conclude that there is any systematic change in Tq or Gf- Stating our results as a 2a 
bound, we have 

<nr 9 yr _1 - (9) 

The main source of error was the fitting errors associated with the determination of 
Nevertheless, the results for the variation of the decay rate of the 56 Co isotope seem 
to be roughly consistent across the five different wavelength bands. 

Most of the statistical pull toward a negative value of comes from two sources, 
SN1995ak and SN1994M, at comparatively high redshifts. The V-band light curves for 
SN 1995ak and the more typical SN 1998aq are shown in figure [3j The two sources at 
the highest z values show markedly more rapid light curve decays than do the lower- z 
supernovae. When sources such as these are used for cosmic distance scale measurements, 
the more rapid decay is removed by applying a stretch factor s to the time axis [2Tj . 
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However, the use of such an empirically determined factor appears to be inappropriate 
here; it could rescale away any actual signal from variations in Gp that was seen from the 
data. If Gp is not varying, some other systematic difference between these sources and 
the others is probably responsible, but the nature of this difference is unknown. 

The second term in eq. (j7|), which comes from the k dependence of ttiq, is considerably 
smaller than the first term, which accounts for its direct dependence on z. Thus, the poorly 
understood physics that causes the light curve to depend on wavelength contributes only 
a small adjustment to the inferred value of -^r-. This was found to be true in all the 
bands. Since the leading order correction due to the /c-dependence is small, as is z for the 
sources we have considered, we do not expect there would be any important changes to ^ 
if higher order in z corrections were taken into account. However, if sources with greater 
redshifts were to be analyzed, the higher order corrections would need to be considered 
and their contributions could be important. 

More sources should be analyzed for this method to produce its best possible result for 
the time variation of the Fermi constant. Although the CfA Supernova Archive contained 
plenty of information for small redshift supernovae (z < 2.4 x 10~ 2 ), more large redshift 
supernovae should also be analyzed in order to obtain better results. At present, the 
10~ 9 yr _1 bound on the fractional time variation can be seen as ruling out roughly a 10% 
change in over about 100 million years of time. Although the late time behavior of 
light curves is more difficult to observe for more distant supernovae, using older sources 
could significantly improve these results. It would also be desirable to have more data 
with well-defined locations for their maxima; light curves with more data points in the 
linear region of interest would also be particularly valuable. 

We found no significant evidence for a change in Gf over the cosmological interval we 
studied, 3.9 x 10~ 3 < z < 2.4 x 10~ 2 . However, if we did see a systematic shift in the 
decay rate of the supernova light curves, it could probably be more reasonably explained 
by something other than a change in the strength of the weak interaction. We have not 
measured Gf (or the half-life of 56 Co) directly. Changes in supernova composition over 
the lifetime of the sample could potentially have affected the decay rate of the light curve, 
changing the way the 56 Co decay products heat the expanding remnant. Statistically 
significant evidence for systematic changes in the supernova light curves over the lifetime 
of the universe would be suggestive, but it would by no means be compelling evidence for 
changing fundamental constants. 

This analysis has determined the rate of change of the dimensional quantity Gf- In 
general, only the time dependence of dimensionless numbers can be defined unambigu- 
ously. A quantity with dimensions must be compared with some similarly dimensioned 
reference unit, which may itself depend on time in a correlated fashion; what reference is 
chosen will affect how the measured quantity appears to depend on time. However, there 
is essentially no problem in this instance. We can calculate how the dimensional param- 
eter T (or Gf) depends on time with no meaningful ambiguity — just as one can observe 
the growth in a child's weight over time without confronting any ambiguity — since any 
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reference with which the dimensional observable is implicitly compared is independent of 
Gp. All possibly relevant time standards with which T could be compared have the same 
dependence on Gp — none. 

In order to extend these bounds to data sets with older, more distant supernovae, 
we would need to account for several additional complications. For very distant sources 
(and especially for late periods along the light curve), the quality of the magnitude data 
diminishes. In addition to limiting the amount of usable data, this could introduce unex- 
pected biases into the data. (For example, light curves with more rapid 56 Co decay rates 
might be underrepresented in the set of usable data, because their luminosity falls below 
the threshold for detection more quickly, providing too few data points to map out the 
relevant region of their light curves accurately.) A more systematic K correction would 
probably also be required at high redshifts, and it would be desirable to make use of more 
detailed models of type la supernovae structure, so as to get at the 56 Co lifetime more 
directly. 

These are not the first published bounds on G F . Although the Fermi coupling con- 
stant cannot straightforwardly be observed through its effects on atomic spectra, it does 
effect nuclear energy levels in a small way. Aside from astrophysical observations of an- 
cient events and precision laboratory experiments, the most successful method for placing 
bounds on the time variations of fundamental constants has been to look at the isotopes 
produced by the natural reactors that operated near Oklo, Gabon about 2 x 10 9 years 
ago. A bound on the fraction change in the dimensionless quantity (3 = Gfw? v at the 
1CT 11 yr^ 1 level was based on the production of samarium isotopes at Oklo [22]. However, 
bounds based on the Oklo data require a very detailed understanding of the operation 
of the reactor, as well as complicated nuclear structure calculations, whose reliability as 
regards the effects of Gp is somewhat questionable. Therefore, our technique, which can 
be extended to older supernovae, may offer the best possibility for placing robust bounds 
on the time variation of the Fermi constant. 
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